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Abstract 

In the last decade Cherenkov telescopes on the ground and space-based y-ray instruments have identified a new class 
of high mass X-ray binaries (HMXB), whose emission is dominated by y rays. To date only five of these systems are 
known. All of them are detected by Cherenkov telescopes in the TeV energy range, while at GeV energies there is 
still one (HESS J0632+057) that has no reported detection with the Fermz'-LAT. A deep search for y-ray emission of 
HESS J0632+057 has been performed using more than 3.5 years of Fermz'-LAT data. We discuss the results of this 
search and compare it to other y-ray binary systems. 
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1. Introduction 

The High Mass X-ray Binaries (HMXBs) are rela- 
tively young (< 10 8 year) Galactic sources composed 
of a massive OB or Be type star and a compact object, 
either a neutron star or a black hole. These systems are 
bright X-rays emitters, generally (but not always) due to 
matter from the massive star accreting onto the compact 
object. 

In the last decade a handful of HMXBs have been 
detected at high (HE; 0.1-100 GeV) or very high- 
energies (VHE; >100 GeV). They are LS I +61°303 
|H11[7]|2|, LS5039 lfTT1l3l. PSR B 1259-63 (ME Ell, 
1FGL J1018.6-5856 |9), Cyg X-3 @, and Cyg X-l 
lfT5ll4D . With the exception of the last two sources, for 
the others the power of their emission is dominant at 
high energies, rather than at longer wave-lengths. For 
this reason they are some times mentioned with the 
name of y-ray HMXBs. Recently a new source has 
been claimed to belong to this class of sources. It is 
HESS J0632+057. In this proceeding we summarize the 
path covered from the discovery of HESS J0632+057 as 
unidentified source, to the claiming of its binarity. Fi- 
nally, we describe the deep search for its emission in 



the GeV energy range, and we discuss the analogies and 
the differences that this particular source shows with re- 
spects the others y-ray HMXBs. 

2. Discovery, and Cherenkov monitoring 

HESS J0632+057 was discovered at VHE energies 
by the High Energy Spectroscopic System (H.E.S.S.) on 
2007 [ 1 3 1 . It is in the Monoceros Loop region, very 
close to the Rosette Nebulae. It was classified as a 
point-like unidentified source. It is interesting to note 
that all the VHE Galactic sources are extended, with 
the exception of the HMXBs. Therefore, the point- 
like feature of HESS J0632+057 was a strong indi- 
cation of its binary nature. Further more, of many 
tens of unidentified sources discovered in the Galac- 
tic Plane H.E.S.S. survey, HESS J0632+057 is one of 
the only two unidentified VHE sources which appears 
to be point-like within the experimental resolution of 
the instrument. The other (HESS J1745-290) is coin- 
cident with the gravitational center of the Milky Way 
02). When HESS J0632+057 was discovered, sev- 
eral plausible association with sources in other wave- 
lengths were discussed. A posteriori, the most relevant 
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Figure 1: Left: Figure from f| 1 91 . X-ray light curve of HESS J0632+057 folded over the period of 321 days. The lower panel shows the hardness 
ratio (2.010.0 keV)/(0.32.0 keV), folded over the same period, and binned at 25 day intervals. Right: Figure from [32). Integral y-ray fluxes above 
1 TeV of HESS J0632+057 from VERITAS (filled markers) and Swift XRT X-ray measurements (0.3-10 keV) in gray. The colors indicate the 
different observing periods of VERITAS. 



was the association with a massive emission-line star 
of spectral type BOpe (MWC 148, see section Rj>. The 
flux and the spectral index of this source measured by 
H.E.S.S. are F(E > ITeV) = 6.4 + 1.5 x 1(T 13 cmT 2 s _1 , 
T = 2.53 + 0.26, respectively. 

Subsiquently to the discovery, VERITAS and 
MAGIC started to monitor HESS J0632+057. It was 
observed by VERITAS for 31 hr in 2006, 2008, and 
2009. The first observation was pointed toward the cen- 
ter of the Monoceros region, (at an angular distance 
of ~ 0°.5 from the HESS J0632+057). The other ob- 
servations were targeted around the reported position 
of HESS J0632+057. During these observations, no 
significant signal in at energies above 1 TeV was de- 
tected, and a flux upper limit of 4.2 x 10~ 13 cirT 2 s _1 at 
99% confidence level |8|. The non detection by VER- 
ITAS during these observations provides variability in 
the VHE flux of HESS J0632+057. This was in favor 
with the binary interpretation of this source. 

Afterwards, the source was detected by MAGIC and 
VERITAS in February 2011, when TeV observation 
were triggered by an outburst in X-rays Il24l [3T1 [TBI . 
The spectrum measured by MAGIC is compatible with 
the one measured by H.E.S.S. Further observations with 
VERITAS have been performed until February 2012 
||32l . In the right panel of figure [TJ are plotted the in- 
tegral y-ray fluxes above 1 TeV of HESS J0632+057 
measured by VERITAS in the different observations. 

3. X-ray follow up, and orbital period detection 

Since its discovery, HESS J0632+057 was followed 
by several X-ray observatories. As first, an observation 



taken on the 2007 September 17 with the EPIC cam- 
era of XMM-Newton lead to the detection of a new X- 
ray source (XMMU J063259.3+054801) within the er- 
ror box of HESS J0632+057, and coincident with the 
massive star MWC 148. This source exhibits a hard 
spectrum, consistent with an absorbed power law with 
r = 1.26 ± 0.04, and the unabsorbed 1 - 10 keV flux is 
(5.3 ± 0.4) x 10~ 13 erg cirT 2 s _1 . 

A long monitoring of this source was planned with 
the Swift X-ray telescope. A first set of data covering 
108 days was analyzed by [23]. They measured flux in- 
creases by a factor ~ 3., and found that X-ray variabil- 
ity was present on multiple timescales including days to 
months. 

The X-ray outburst on February 2011 that was an- 
nounced by [24 1 triggered the observation of VERITAS 
and MAGIC (see section [2]), as well as an observation 
with Chandra requested by [40|. Comparing the re- 
sults of the Chandra data analysis with those of the 
first observation of XMM-Newton [40 1 found also a a 
spectral variability in addition to the flux one. Indeed, 
the spectral index measured was 1.61 + 0.03. The flux 
and spectral variability are both characteristic features 
of the HMXBs, but in contrast with what observed for 
the other TeV binaries, in this source the higher the flux 
the softer the X-ray spectrum. 

Finally, the detection of the orbital periodicity of 
321+5 days was claimed by |fT9l . The periodical X- 
ray flux modulation was the result of the analysis of the 
long set of Swift observations of the source from 2009 
January to 2011 March. In figure [TJ (left plot) is shown 
the folded X-ray light curve of HESS J0632+057. It 
shows a narrow high peak at phase 0.3 followed by a 
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dip feature, and centered at phase 0.75 there is a mod- 
erate broad peak. The XMM-Newton observation fall in 
the dip feature, where the spectrum is harder, while the 
Chandra observation was taken during the high peak. 
Phase 0.0 in this plot is arbitrary set at the beginning of 
the Swift observations. 



4. Optical and Radio counterparts 



MWC 148 is the companion massive star in the bi- 
nary system HESS J0632+056. It is a Be star, as in 
the case of the HMXB systems LS I +6T303, and PSR 
B 1259-63. Before the detection of the orbital period 
in X-rays, IfTTI analyzed optical spectra of the star ac- 
quired to investigate the stellar parameters, the proper- 
ties of the Be star disk, and evidence of the binarity. 
They derived a mass of 13.2-19.0 M , and a radius of 
6.0-9.0 R . Fitting the spectral energy distribution they 
found a distance between 1.1 and 1.7 kpc. Since their 
dataset was too short (~ 1 month) it was not possible 
to detect the orbital period of the system. On the other 
hand, [21 1 collected optical spectra observations span- 
ning for 4 years from 2008 October to 201 1 May. They 
calculated the orbital parameters of the binary systems, 
finding that it is highly eccentric with e = 0.83 ± 0.08. 
Furthermore, the orbital solution implies that the high 
peak in the X-ray light curve is delayed ~ 0.3 orbital 
phases after the periastron passage, similar to the case 
ofLSI+61°303. 

In the radio regime [43 1 detected with the Giant Me- 
trewave Radio Telescope (GMRT) and the Very Large 
Array (VLA) observations at 1.28 and 5 GHz, respec- 
tively, a point-like, variable radio source at the posi- 
tion of MCW 148. 01 observed HESS J0632+057 
with the European VLBI Network (EVN) at 1.6 GHz 
in two epochs: during the January/February 2011 X- 
ray outburst and 30 days later. In the first epoch the 
source appears point-like, whereas in the second one it 
appears extended with a projected size of ~75 AU as- 
suming a distance to the system of 1.5 kpc. The bright- 
ness temperature of 2xl0 6 K at 1.6 GHz together with 
the negative spectral index around -0.6 point to non- 
thermal synchrotron radiation producing the radio emis- 
sion 11431 [34l . The detected projected displacement of 
the peak of the emission of 21 AU in 30 days, the mor- 
phology and the size of HESS J0632+057 are similar to 
the ones observed in the well-established y-ray binaries 
PSR B1259-63, LS 5039, LS I +61°303. 



Figure 2: Background subtracted count map above 1 GeV of the 10° X 
10° sky region centered on HESS J0632+057. The pixel size is 0.1° X 
0.1°. A Gaussian smoothing with 3 pixel kernel radius is applied. 
The 2FGL sources are labeled with a cross, HESS J0632+057 by the 
diamond, and PSR J0633+0632 by the circle. 

5. Search for GeV emission 

All the y-ray HMXBs have a bright emission in the 
GeV energy range, but HESS J0632+057 is not present 
in the Second ferm/-LAT Source Catalog ll35ll (2FGL 
hereafter). Therefore, we performed a fine analysis with 
the Fermi-LP^L data in order to deeply search for GeV 
emission. A dataset of ~ 3.5 years from 4 th August 
2008 to 3 rd March 2012 has been analyzed, that is al- 
most twice the dataset used for the 2FGL. 

HESS J0632+057 lies at b = -1.44° from the Galac- 
tic plane, in the active region of the Monoceros loop. 
The region at GeV energies appears as shown in Fig- 
ure [2] where the sources from the 2FGL catalog are la- 
beled with a green cross, the yellow one in the center 
marks the position of HESS J0632+057, while in cyan 
is marked the pulsar PSR J0633+0632. It is evident that 
the pulsar emission dominates the region aroung HESS 
J0632+057. Their angular distance is only 0.78°. 

5.7. PSRJ0633+0632 

PSR J0633+0632 is a bright radio-quiet y-ray pulsar 
discovered in the first six months of the Fermi mission 
ID. Indeed, it is listed in the first Fermi-LAT catalog 
of y-ray pulsars [5]. It has a flux of 8.4 + 1.2 x 10~ 8 
ph cirT 2 s _1 , and the pulse profile has the characteris- 
tic two narrow peaks separated ~ 0.5 in phase. The 
weighted pulse profile of PSR J0633+0632 is shown in 
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Figure3: Weighted pulse profile of PSR J0633+0632 with- 3.5 years 
collected Femi-LAT data. The weight associated to each event count 
corresponds to the probability that the event is emitted by the pulsar, 
rather than by a nearby source or by the diffuse emission. The off- 
pulse, and the intra-peaks phases are defined by mean of a Bayesian 
block algorithm (green line): (0.66 - 0.98) and (0.21 - 0.49), respec- 
tively. 



Figure 4: Diffuse subtracted TS map of the 5° x5° sky region centered 
on HESS J0632+057. The maps are calculated for E > 300 MeV. 
The 2FGL sources are labeled with a cross, HESS J0632+057 by the 
yellow circle. 



5.2. Upper limits of HESS J0632+057 



figure [3] A weight was associated to each event, corre- 
sponding to the probability that the event is emitted by 
the pulsar, rather than by a nearby source or by the dif- 
fuse emission. We defined the off-pulse phases of PSR 
J0633+0632 (0.66 - 0.98) applying a Bayesian block al- 
gorithm adapted for weighted-counts light curves [20]. 
In figure[3]the green line shows the calculated Bayesian 
blocks. The off-pulse is defined as the lowest block 
reduced by 10% of its total length at each edge. As 
well, we defined the intra-peaks or bridge phases (0.21 
- 0.49). 

In order to avoid the contamination from the strong 
emission of PSR J0633+0632, we carried on the analy- 
sis of HESS J0632+057 selecting the off-pulse and the 
intra-peaks phases of PSR J0633+0632. 

With this phase selection the regioin around HESS 
J0632+057 appears quite crowded, and not perfectly 
modeled by the 2FGL sources. This is mainly due to 
the fact that the 2FGL catalog is built with 2 years of 
Fermz'-LAT data, while in this analysis the dataset is al- 
most doubled (~ 3.5 years). 

Figure [4] shows the region out of the pulsar peaks. 
The 2FGL sources are labeled with a cross, and HESS 
J0632+057 by the yellow circle. Before to analyze di- 
rectly HESS J0632+057 we as first applied an itera- 
tive method to reach a good modeling of the region 
around it. The method developed is described in detail 
in ll20l . It accounts for the excesses in the map adding 
new point-like sources in the model of the region. 



Once we obtained a good model for the region, we 
came back to focus specifically on the analysis of HESS 
J0632+057. We found no significant detection of the 
GeV emission of HESS J0632+057. Then, we derived 
the 95% flux upper limit for this source. The Helene's 
Bayesian method 11251 was adopted to evaluate them. 
We calculated the flux upper limit for energies E > 100 
MeV of Fioo < 3.0 x 10~ 8 cirT 2 s _1 . 

5.3. Folded analysis and variability 

In the analysis described so far we have considered 
HESS J0632+057 as a steady source, and we concluded 
that its average y-ray emission is below the sensitivity of 
Fermz'-LAT. But it is still possible that the system have 
a brighter emission during a specific part of the orbit, 
like for the case of the HMXB system PSR B 1259-63 
(ED, ED). In order to investigate this hypothesis, a fur- 
ther analysis was performed folding the Fermz'-LAT data 
with the orbital period of 321 days, and setting the phase 
0.0 to the epoch 54857 MJD |[T9). We subdivided the or- 
bit in 8 equally spaced phase bins, and for each of them 
we calculated the significance of the GeV emission. The 
results are showed in figure[5] We found that none of the 
bins have a significacnce above the detection threshold. 
Finally, we also searched for flares or high flux states 
of HESS J0632+057, calculating light curves with dif- 
ferent time binning scales, from days to months. For 
none of the time scales investigated we got a significant 
signal. 
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Figure 5: The TS light curve of HESS J0632+057 calculated folding 
the Fermi-LAT data with the orbital period of 321 days. The TS values 
of all the phase bins are below the detection threshold (TS= 25). 



6. Conclusions 

In this proceeding we have summarized all the steps 
leading to the claim that HESS J0632+057 is a y-ray 
HMXB system. Fermi-LAT data has been analyzed 
with the aim to search for the emission in the GeV 
energy range of HESS J0632+057. The analysis per- 
formed needed particular attention especially for the 
presence of a very nearby y-ray pulsar PSR J0633+0632 
(section 5.1). As first, we selected the events out of 
the peaks of the pulsars. Then, the region around the 
source was modeled adding new point-like sources in 
correspondence of unmodeled excesses with an itera- 
tive method. Despite the good modeling of the region, 
we did not detected any significant signal of y-ray emis- 
sion from HESS J0632+057. We calculated the upper 
limit for its permanent GeV emission, Fioo < 3.0 x 10~ 8 
cirT 2 s . We also searched for GeV emission from a 
restricted portion of the orbit phases, as well as for plau- 
sible flares, or changes of state in its flux. 

The lack of significant signal from all these searches 
is difficult to understand if compared with the behav- 
ior of the other y-ray HMXBs in the same energy 
range. Similar to HESS J0632+057, the massive 
stars in LS I +61° 303 and PSR B 1259-63 systems 
are Be star, that are characterized by the presence of 
a dense circumstellar disk. The distances of these 
three HMXBs are comparable: ~ 1.5 kpc for HESS 
J0632+057, ~ 1.9 kpc for LS I +61°303, and ~ 2.3 
kpc for PSR B 1259-63. In contrast their orbital 
periods are significantly different. LS I +61° 303 has 
the shortest orbital period of about one month (26.5 
days). As consequence, the systems is also the most 
compact among the three. In the GeV energy range, 
LS I +61°303 has a continuous emission modulated by 
its orbital motion. Most probably the GeV emission 
is present in every phase of the orbit because it is so 
small that the compact object is always interacting 



with the circumstellar disk, or with the stellar wind. A 
different case is that of PSR B 1259-63 with its very 
large orbital period of about 3.5 years. The orbit of this 
system is highly eccentric, with eccentricity e ~ 0.87. 
This makes that when the compact object is close to 
the periastron, it pass through the circumstellar disk. 
Only during this passage PSR B 1259-63 emits GeV 
and TeV y-rays. HESS J0632+057 has an orbital 
period of about 1 year, that is in between those of 
LS I +61°303 and PSR B 1259+63. The eccentricity of 
the orbit of HESS J0632+057 is also very high, with 
e ~ 0.83. We expected a GeV emission at least during 
the periastron passage, but we did not detected it. This 
can be attributed to different factors. An important role 
can be played by the features of the circumstellar disk, 
like its size, its density, and its inclination respect to 
both the observer line of sight and the plane of the orbit. 
On the other hand, the differences among LS I +61° 303 
PSR B 1259+63, and HESS J0632+057, can also be 
attributed to the features of the compact objects, like 
the intensity of their relativistic winds. Or they could 
be caused by the different nature of the compact object 
itself, either a pulsars (as for PSR B 1259+63) or a 
black hole. All these hypothesis are currently under 
investigation. 
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